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Abstract: 

We report on a search for mutually unbiased bases (MUBs) in 6 dimensions. 
We find only triplets of MUBs, and thus do not come close to the theoretical 
upper bound 7. However, we point out that the natural habitat for sets of 
MUBs is the set of all complex Hadamard matrices of the given order, and we 
introduce a natural notion of distance between bases in Hilbert space. This 
allows us to draw a detailed map of where in the landscape the MUB triplets 
are situated. We use available tools, such as the theory of the discrete Fourier 
transform, to organise our results. Finally we present some evidence for the 
conjecture that there exists a four dimensional family of complex Hadamard 
matrices of order 6. If this conjecture is true the landscape in which one may 
search for MUBs is much larger than previously thought. 
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1. Introduction 

This is a paper on 6 by 6 matrices. Its scope would therefore seem rather 
limited, but we will address two much discussed open problems. The first 
concerns the classification of complex Hadamard matrices, and the second 
the existence (or not) of complete sets of mutually unbiased bases. These 
problems are connected to each other, they are of considerable interest in 
quantum information theory, and they have a long history. The study of 
complex Hadamard matrices was begun by Sylvester [1], and by Hadamard 
himself [2], while the second problem has appeared (in various guises) in 
quantum theory [3, 4], operator algebra theory [5], Lie algebra theory [6], 
and elsewhere [7, 8]. It would seem as if all questions concerning single digit 
dimensions should have been settled by now, but this is not the case. 

A complex Hadamard matrix is, by definition, a unitary matrix all of 
whose matrix elements have equal modulus. Much of the mathematical lit- 
erature concerns the special case of real Hadamard matrices, but from now 
on we will take "Hadamard matrix" to refer to a complex Hadamard matrix. 
An obvious question concerns the classification of all N x N Hadamard ma- 
trices. It has been settled for N < 5 [9], but for all larger values of N it is 
open [10, 11]. 

The columns of a unitary matrix define an orthonormal basis in Hilbert 
space. The basis defined by an Hadamard matrix has the property that 
the modulus of the scalar product between any one of its vectors and any 
vector in the standard basis (whose vectors has only one non-zero entry) is 
always equal to 1/ y/N. Again by definition, we say that these two bases are 
mutually unbiased, or MUB for short. One can now ask how many bases, 
with all pairs being MUB with respect to each other, one can find in a Hilbert 
space of a given dimension N. It is known that there can be at most N + 1 
MUBs, and when N is the power of a prime number the answer is N + 1 [4]. 
For other values of iV the question is open. For N = 6, all that is known is 
that the number of MUBs is at least 3 and at most 7 [12]. 

The set of Hadamard matrices provides a landscape where sets of MUBs 
live. In fact classifying the set of Hadamard matrices is equivalent to clas- 
sifying the set of ordered MUB pairs, up to natural equivalences [13]. Once 
such a classification has been carried out we proceed to list all MUB triplets, 
all MUB quartets (if any), and so on. However, let us admit at the outset 
that we have carried through this strategy only piecemeal — and for N = 6 
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we did not find any MUB quartets. 

In section 2 of this paper we describe the known Hadamard matrices 
of order 6, and sort out some problems concerning equivalences between 
them. In section 3 we recall some facts about MUBs, and decide when 
sets of MUBs should be regarded as equivalent. In section 4 we define a 
useful notion of distance between orthonormal bases. The distance attains 
its maximum when the bases are unbiased. In section 5 we describe a search 
for all MUBs whose vectors are composed of rational roots of unity of some 
modest orders; the most interesting of our choices is the 24th root. In section 
6 we describe the set of all bases that are MUB with respect to both the 
standard and the Fourier bases [14], and show how this set can be understood 
from properties of the discrete Fourier transform. In section 7 we describe 
MUBs of a related type, where the Hadamard matrices have a particular 
block structure. Finally, in section 8 we give some arguments suggesting that 
there exists a four dimensional family of Hadamard matrices — the largest 
known family has two dimensions only. Readers interested in Hadamard 
matrices only can skip sections 3-7, but our main contention is that readers 
interested in MUBs should skip nothing. Section 9 states some conclusions. 

Notation: z is a complex number, and z its complex conjugate. The ratio- 
nal root of unity e 2m ^ n is denoted u. Its integer powers are called nth roots. 
A special case is q = e im l N ; where N is the dimension of the Hilbert space. 
A matrix element of the matrix M is denoted M ao , M T is the transposed 
matrix, and is the adjoint. The index a usually runs from to 5. 

2. Complex Hadamard matrices 

According to our definition an Hadamard matrix H is a unitary matrix whose 
matrix elements obey 



An example that exists for all N is the Fourier matrix F, whose matrix 
elements are 




< a,b < N - 1 . 



(1) 



1 



,ab 




(2) 



F a b 




q 
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Two Hadamard matrices Hi and if 2 are called equivalent [9], written H\ m 
H2, if there exist diagonal unitary matrices D\ and D 2 and permutation 
matrices Pi and P 2 such that 

#1 = D1P1 • H 2 ■ P 2 D 2 . (3) 

That is to say, we are allowed to rephase and permute rows as well as columns. 
This is an equivalence relation, and we are interested in classifying all the 
equivalence classes for a given matrix size N. 

We will usually present our Hadamard matrices in dephased form, which 
means that all elements in the first row and the first column are real and 
positive. This can always be achieved using diagonal unitaries, but it does 
not fix the equivalence class completely. If the Hadamard matrix is written 
in any other form it is said to be enphased. 

For N = 2, 3, and 5, the Fourier matrix is unique, in the sense that it 
represents the only equivalence class of Hadamard matrices [9]. For N = 
4 there exists a one parameter family of equivalence classes, including the 
Fourier matrix as well as a real Hadamard matrix [2]. This family exhausts 
the set of equivalence classes when N = 4. Let us remark that continuous 
families appear also when N — 7 [11], so that their existence does not hinge 
on N not being a prime number. 

We will print concrete Hadamards matrices in boldface. When N = 6 the 
following representatives of different equivalence classes are known to us: 

• A two parameter family F(xi, x 2 ), including the Fourier matrix F(0, 0). 

• The transpose F T (xi,x 2 ) of the above. 

• A circulant matrix C found by Bjorck [16], and its complex conjugate. 

• A one parameter family D(x), including a matrix D composed of fourth 
roots of unity, called Diva's matrix [10]. 

• A one parameter family B(#) that interpolates between C and D [17]. 

• Tao's matrix S, composed of third roots of unity [18, 19]. 

Attribution may be difficult; the Dfya family, or representatives thereof, was 
discovered several times [20, 9, 12]. With one exception, namely the family 
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B(0), the known continuous families are of the special kind called affine 
families [11], which means that some of the matrix elements can be multiplied 
with free phase factors in such a way that the matrix remains Hadamard. 
We now go through the items in our list. The Fourier family is explicitly 



F(x 1 ,x 2 ) 



qzi 
q 3 z 1 



q 5 zi 



1 1 

2 S 

q z 2 q 



Z2 

q 2 

4 

q z 2 



q' 



1 

4 

q zi 

q 2 

Zl 

q 4 

q 2 zi 



q 5 z 2 

q 4 

q 3 z 2 

q 2 

qz 2 



zi = e 27riXl 



z 2 = e 27TiX2 



(4) 



The two free parameters arise essentially because a six dimensional space can 
be written as a tensor product; see section 7. This is an affine family, since 
there are no restrictions on the phase factors z\ and z 2 . 



Figure 1: The affine family F(xi,a?2) is divided into 144 equivalent triangles of 
equal area, as described in the text. The lattice shown in the background are the 
6th roots. 

There are equivalence relations connecting different matrices within this 
family. They form a discrete group, and thus they will affect the size but not 
the dimensionality of the affine family. Let us think of the parameter space 
as a square in the (xi,x 2 ) plane (see fig. 1). By inspection one finds cyclic 
translation groups of order 6: 
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F(x 1 ,x 2 ) « F(xi + 2/6, x 2 + 1/6) w F(xi + 4/6, x 2 + 2/6) « F(x 1? -rr 2 ), 

(5) 

F(xi, ar 2 ) « F(xi + 1/6, x 2 + 2/6) w F(xi + 2/6, x 2 + 4/6) w F(-xi, x 2 ). 

They divide the square into 12 copies of a fundamental region that is a paral- 
lellogram. There are further equivalences corresponding to transformations 
of order 2. By inspection 

F(x 1 ,x 2 ) « F(x 2 ,xi) . (6) 

If we permute the rows and then dephase the row that ends up in the first 
place, we find that 

F(x u x 2 ) w F(-x u -x 2 ) ■ (7) 

As a result, our parallellogram will be divided into four equivalent triangles. 
Finally we can permute the columns and dephase the column that ends up 
in the first place. This gives a final cyclic equivalence group of order 3: 

F(x!,x 2 ) F(x 2 - xi, -xi) ~ F(-x 2 , xi - x 2 ) (8) 

Considered as transformations of our square, these last transformations are 
not isometries but they do preserve area. Hence the triangles that we had in 
the previous step are themselves divided into three triangles of equal area, one 
of which is our final fundamental region of inequivalent Hadamard matrices. 
A convenient choice of fundamental region is a triangle with corners at (0, 0), 
(1/6,0) and (1/6, 1/12); the original square is covered by 144 copies of this 
region. The family F T (xi,x 2 ) works (at this stage) in an entirely analogous 
way. 

Bjdrck's circulant Hadamard matrix is 

1 id —d —i —d id 

id 1 id —d —i —d 

£, _ — d id 1 id —d —i 

—i —d id 1 id —d 

—d —i —d id 1 id 

id —d —i —d id 1 



dd = 1 . (9) 
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See section 6 for further discussion of circulant Hadamard matrices. The 
complex number d has modulus unity and is 



d = ^^ + t\jY- d 2 -{1-V3)d+1 = . (10) 

No affine family stems from this matrix [11]. We will have more to say about 
this matrix later; for now let us just mention that all Hadamard matrices 
equivalent to a circulant matrix have been listed by Bjorck and coauthors 
[16], for N < 8. 

The Difa family is 
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z = e 2vix 



(11) 



There are equivalence relations within this family too. In fact 

D(x) « T>(x + 1/2) « D(-x + 1/4) . (12) 

Hence we can take — 1/8 < a: < 1/8 without loss of generality. The matrices 
D(x) can be written on block circulant form; see section 7. 

The Hermitian family B(0) was found very recently, by Beauchamp and 
Nicoara [17], and consists — up to equivalences — of all Hermitian Hadamard 
matrices of order 6. Explicitly it is 
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where (x, y, z, t) are complex numbers of modulus one, related by 

y = e %e t = xyz 



(13) 



(14) 



7 



z - 1 + 2y ~ y2 (15) 
l + 2y + y 2 ±V2^/l + 2y + 2y* + y4 

x = TT^? ' (16) 

with 9 providing a free parameter. The two branches of the square root lead 
to equivalent families. This is an example of a non-affine family (and the 
only such example known to us for N — 6). 

The phase 9 cannot be chosen arbitrarily; an interval around 9 = is 
excluded. The Hadamard property requires that 

\y\ 2 = l =>■ \x\ 2 = \z\ 2 = \t\ 2 = 1 . (17) 

On examining this point one sees that this restricts the allowed values of 
y = cos 9 + isin#. There are no difficulties with z, but the numerator of x 
becomes 



l + 2y + y 2 ±y/2yjl + 2y + 2y 3 + y 4 = 2y(l + cos0±V2 cosfl + 2 cos 2 6 - 1). 

(18) 

The absolute value of this expression hinges crucially on the sign inside the 
square root; the absolute value of x will be equal to one if that sign is negative. 
This restricts the phase to 

cos 0<^Lli. (19) 

The allowed range of 9 has end points, with vanishing square root, at y — 
—d and y = —d. The curve is not closed; its end points correspond to 
permutations of Bjorck's matrix C and its complex conjugate. This is also 
the case if y = d 2 or y = d 2 . If y = — 1 or y = ±i we obtain permutations of 
Dita's matrix. 

Tao 's matrix never appeared in our searches for MUBs, so we do not give 
it explicitly here. Since Tao's matrix is composed of 3d roots only, one may 
ask whether there are any N = 6 Hadamard matrices based only on 5th or 
7th roots. Using a computer program to be described in section 5 we have 
shown that no such Hadamard matrix exists. 



8 



3. Preliminaries on MUBs 

We now turn to Mutually Unbiased Bases, or MUBs. First we recall some 
definitions. Fix an orthonormal basis \e a ) in an N dimensional Hilbert space. 
A unit vector |/) is said to be unbiased with respect to the fixed basis if for 
all a 

\(e a \f}\ 2 = ^. (20) 

The important thing is that the right hand side is constant; its value is a 
consequence. A pair of orthonormal bases is said to be an unbiased pair if all 
the vectors in one of the bases are unbiased with respect to the other. One 
can go on to define larger sets of mutually unbiased bases, or MUBs, in the 
obvious way. It is known that the number of MUBs one can find is bounded 
from above by N + 1, but it is not known if this bound can be achieved unless 
N is a power of a prime number. A set of N + 1 MUBs, if it exists, is known 
as a complete set. 

If one member of a set of MUBs is represented by the columns of the 
unit matrix, all the other bases must be represented by the columns of a 
set of Hadamard matrices. A complete set of N + 1 MUBs exists if we can 
find N enphased Hadamard matrices representing bases that are MUB with 
respect to each other. Altogether then we have N Hadamard matrices that 
are said to be Mutually Unbiased Hadamards, or MUHs [11]. Note that two 
Hadamard matrices Hi , H 2 are unbiased if and only if 

H\H 2 = H 3 , (21) 

where H 3 is an Hadamard matrix too. We say that a set of MUHs all of whose 
members are equivalent in the sense of section 2 is homogeneous, otherwise 
it is heterogeneous. The standard construction of complete sets of MUBs in 
prime power dimensions [4] gives a homogeneous set of N MUHs. 

When are two pairs of MUBs equivalent to each other? Let (M , Mi) 
denote an ordered pair of MUBs, and {M , Mi} an unordered pair, with each 
basis represented as the columns of a unitary matrix. We are not interested 
in the order of the constituent vectors, nor in the phase factors multiplying 
these vectors. Hence, with P a permutation and D a diagonal unitary matrix, 

(M PD, MiP'D') w (M , Mi) . (22) 
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We also declare that an overall unitary transformation of the Hilbert space 
is irrelevant, so that 

(C/M , UM ± ) « (M , Mi) . (23) 

Then the pair can always be represented in the form (1, H), where H is 
an Hadamard matrix. It follows that the equivalence relation (3), used for 
Hadamard matrices, is the correct equivalence relation also for ordered MUB 
pairs; 

(1, H ± ) « (1, H 2 ) <£> H x kH 2 . (24) 

From eq. (23) we know that (1, H) ps (H\ t). For unordered MUB pairs this 
means that 

{1,H} . (25) 

Therefore, the condition for unordered pairs to be equivalent is 

either Hi ~ H 2 

{l,Hi}^{l,H 2 } & (26) 

or Hi ps H\ 

The discussion can be extended to larger sets of MUBs, ordered or unordered 
—although the equivalences will become harder to check as the number of 
MUBs increases. 

For the Hadamard matrices we know, it is true that 

[F(xi,x 2 )]^F T (xi,x 2 ) (27) 
[D(x)] f w D(-x) (28) 

C 1 " w C S f « S , (29) 

and finally the Beauchamp-Nicoara family is Hermitian by construction. 
Thus the set of unordered MUB pairs is smaller than the set of inequiv- 
alent Hadamard matrices. In particular, {1, F(xi, x 2 )} and {t,F T (xi,x 2 )} 
are equivalent when considered as unordered pairs. 
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Let us return to the MUBs themselves. It is useful to think of them as 
sets of density matrices, rather than as sets of vectors in Hilbert space. Recall 
that a unit vector |e) in an N dimensional Hilbert space corresponds to a 
projector |e)(e|, which is an Hermitian matrix of trace unity, and as such can 
be regarded as a vector in an N 2 — 1 real dimensional vector space whose 
elements are matrices. The origin of this vector space is naturally chosen 
to sit at the matrix = j^l, and then its vectors are traceless matrices. 
Explicitly the correspondence is 



|e> - e = y^(|e)< e |-A0 . (30) 
The vector e is a unit vector with respect to the scalar product 

e • f = ^Tr ef . (31) 
The distance between two matrices A and B is given by 

d 2 (A,B) = ^Tr(A-B) 2 . (32) 

This is the Euclidean Hilbert-Schmidt distance. 

It is important to realize that although any unit vector in Hilbert space 
gives rise to a unit vector in the larger space, it is only a small subset of the 
latter that can be realized in this way. The convex cover of this small subset 
is the set of all density matrices (or quantum states). An orthonormal basis 
in the Hilbert space corresponds to N vectors that form a regular simplex 
in the larger space. This simplex spans an N — 1 dimensional plane through 
the origin. It is moreover easy to see that 

l(e|/>| 2 = ^ * e-f = 0. (33) 

Therefore, if the unit vectors belong to a pair of MUBs, the corresponding 
vectors in the large vector space are orthogonal. It follows that a pair of 
MUBs span two totally orthogonal (N — l)-planes through the origin. This 
fact is at the bottom of the reason why MUBs are useful in quantum state 
tomography [4]. And we see immediately that there can be at most N + 1 
MUBs, since there is no room for more than N + l totally orthogonal (TV — 1)- 
planes in an N 2 — 1 dimensional space. 
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4. Interlude: a distance between bases 



We need a distance between bases in Hilbert space, such that it becomes 
maximal if the bases are MUB. It should be natural and easy to compute, 
but we do not insist on any precise operational meaning for it. 

As explained in section 3, a basis in an N dimensional Hilbert space spans 
an (N — l)-plane in a real vector space of N 2 — 1 dimensions. There is a 
standard way to define a distance between such planes, namely to regard 
them as points in a Grassmannian, in itself embedded in the surface of a 
sphere within an Euclidean space of a quite high dimension [21]. The proce- 
dure is analogous to the one we followed when we transformed our Hilbert 
space vectors into points in a space embedded within the space of traceless 
matrices. The details are as follows. Starting from a basis \e a ) in Hilbert 
space, use the recipe given in eq. (30) to form the N vectors e a . Expand 
these vectors relative to a basis. Then form the (N 2 — 1) x N matrix 



N - 1 

5 = y n ^[e 1 e 2 ... e N }. (34) 

It has rank N — 1. Next introduce an (N 2 — 1) x (iV 2 — 1) matrix of trace 
N — 1, projecting onto the (N — 1) dimensional plane spanned by the e a : 



r, T N ~ 1 r 

P = B B = — — [ei ... e N ] 



T 



(35) 



It is easy to check (through acting on e a say) that this really is a projector. 
The chordal Grassmannian distance between two (N — l)-planes is defined 
in terms of the corresponding projectors as 

D " {Pl ' P2) = 2(N-l) Tl{Pl " P2)2 = 1 " iv~T TrPlP2 • (36) 

There is an analogy to how the density matrices were defined in the first 
place, and to the Hilbert-Schmidt distance between them. 

We can think of the projectors P as points on the surface of a sphere in 
R M , where — as it happens — 

N 4 -N 2 -2 
M = . (37) 
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This explains the name "chordal distance". It is however important to re- 
alize that the Grassmannian of (N — l)-planes forms a very small subset of 
this sphere. Its dimension is N(N — l) 2 . And then bases in Hilbert space 
correspond to a small subset of the Grassmannian. 

The chordal distance has been used in studies of packing problems for 
planes and subspaces [22] . It has a number of advantages when compared to 
more sophisticated distances, such as the geodesic distance within the Grass- 
mannian (which would be analogous to the Fubini-Study distance between 
pure quantum states). It is useful to know that the chordal distance can be 
written as a function of the principal angles Oi, namely 

1 N-l i N-l 

Dl = 1 - j^-r E ™s 2 ^ ~E sin 2 ^ • (38) 

JV 1 i=l JV 1 i=l 

Here the first principal angle is defined as 

cos#i = maxui ■ vi , (39) 

where ui and vi are unit vectors belonging to the respective (iV— l)-planes, 
chosen so that their scalar product is maximized. The second principal angle 
is defined by maximizing the scalar product between unit vectors in the 
orthogonal complements to Ui and Vi, and so on. It is then clear that 

O^DlKl-^, (40) 

where k is the dimension of the intersection of the two planes. The distance 
is maximal if and only if the (N — l)-planes are totally orthogonal. 

Now consider two (N — l)-planes spanned by vectors corresponding to 
two bases \e a ) and \f a ) in Hilbert space. Working through the details, one 
finds that the distance squared between the bases is 

1 N-l N-l , 1 s 2 

^^) = 1 -oEENl J -i • ( 41 ) 

JV 1 a=0 6=0 V JV 7 

Thus, between bases in Hilbert space, 

< D 2 C < D 2 ia x = 1 • (42) 
The distance attains its maximum value if and only if the bases are MUB. 
A set of MUBs forms a regular equatorial simplex on the sphere in R M , 



13 



although there will be many regular equatorial simplices that do not arise in 
this way. 

What is D 2 C on the average, for two bases picked at random in Hilbert 
space? To answer this question we represent one basis by the standard basis. 
Then we choose a vector at random according to the Fubini-Study mea- 
sure, choose a vector in its orthogonal complement again according to the 
Fubini-Study measure (in one dimension lower), and so on until we have a 
complete basis. The resulting measure is a measure on the flag manifold 
U(N)/[U(1)] N . In practice the calculation is simple; the average is 



c 



Dl = 1 



1 



N * a b 



rrEEIJ^I 2 -^) > • (43) 




Using the linearity of the average, together with the fact that the N 2 terms 
in the sum must have equal averages, we can rewrite this as 

D -) = 1 -i&i((^ 2 -h) 2 ) • (44) 

Hence it is enough to calculate the average of a function of the modulus of a 
component of a random vector, using the Fubini-Study measure. How to do 
this is described elsewhere [23, 21]; the answer we arrive at is 

Dl) = — — - . 45 

c / N + l v ' 

For N = 6 we have (D 2 ) = 0.86; as iV grows the average distance squared 
approaches the maximum value 1. Note that the average depends smoothly 
on N, whatever the maximal number of MUBs in N dimensions may be. 

One can ask for a function of iV + 1 points on the sphere in R M , whose 
maximum is attained when the points form a regular equatorial simplex. 
Given that N + 1 < M it happens that 

N+l N+l 

/-EE^,^), (46) 
i=l j=l 

is such a function. To see this, use the Euclidean norm on the space in which 
the Grassmannian of (N — l)-planes is embedded, and normalise it so that 
the projectors correspond to N + 1 unit vectors Ej. Then 
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^||E i -E,|| 2 = iV(iV+l)-2^E J -E j . (47) 

i<j i<j 

On the other hand 

iV + 1 + 2^Ei • Ej = II^EiH 2 > . (48) 

i<j i 

A minimum of the last expression corresponds to a maximum of the function 
/. A regular equatorial simplex clearly saturates the inquality (48), given 
that N + 1 is smaller than the dimension of the space that we are in. This 
proves our point. There are other configurations besides the desired one that 
also saturate the bound. Since they do not arise from bases in the underlying 
Hilbert space we can ignore such configurations. Still we do not claim that 
/ is the most useful function of its kind. 

Equipped with our notion of distance (so that we can quantify our fail- 
ures), we can look for complete sets of 7 MUBs in C 6 . The natural way to 
proceed, given the results of this section, is to maximise the function 

/ = £XX 2 (^) = EE fl"^ EE (l(e«l4 J) >| 2 -^) 2 ) , (49) 

i=lj=l i=lj=l\ °a=0fe=0 V ' ) 

with the understanding that the \e$) are orthonormal bases in Hilbert space. 
There is no a priori reason to believe that the upper bound on / can be 
attained, because we are confined to a small subset of all possible (N — 1)- 
planes. Still, it might be interesting to find the "best" solution in this sense. 
A similar but more sophisticated procedure has been successfully used to find 
a special kind of overcomplete bases known as SIC-POVMs [24]. We have 
not attempted such a calculation however. 

5. MUBs composed from rational roots of unity 

In our search for MUBs we rely on the classification of ordered MUB pairs 
through Hadamard matrices. Without loss of generality we assume that the 
first MUB is represented by the standard basis 1, and call it (MUB) . Then 
we choose a representative of an equivalence class of Hadamard matrices, 
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and call it (MUB)i. Assume it to be given in dephased form. All bases that 
are MUB with respect to the first two must then be represented by enphased 
Hadamard matrices, and we proceed to look for them. 

For N < 5 there exist maximal sets of N + 1 MUBs, and these sets are 
unique up to an overall unitary transformation. Given that (MUB)! is not 
unique for N = 4, this is perhaps a little surprising. But for N = 4 it 
is possible to find inequivalent and incomplete sets of MUBs, that cannot 
be extended to complete sets [25]. To get the complete set one must start 
with the real Hadamard matrix, rather than the Fourier matrix, as (MUB)i. 
Given that Hadamard matrices of order 6 are highly non-unique, the question 
where to place (MUB)i becomes non-trivial for N = 6. 

All known complete sets of MUBs are built from vectors all of whose com- 
ponents are iVth or 2iVth roots of unity, depending on whether the dimension 
N is odd or even [4, 26]. We therefore made a program that lists, for N = 6, 
all orthonormal bases whose vectors are composed of 12th roots. They are 
candidates for (MUB)i. Next, for each (MUB)i the program lists the set of 
all vectors composed from 12th roots and unbiased with respect to it. The 
program also lists all orthonormal bases that can be constructed from each 
set of unbiased vectors. They are candidates for (MUB) 2 . Finally we prune 
the list so that it contains only inequivalent Hadamard matrices as (MUB)i, 
and if there are several (MUB)2 we compute the distance between them to 
see if we get any sets of four MUBs in this way. Note that this means that 
the part of the landscape of Hadamard matrices we search in consists of the 
Fourier and Di$a families, the Tao matrix, enphased versions of these, and 
(probably) nothing more. The analogous calculation when N is a power of 
a prime would have found all known complete sets of MUBs, but for N = 6 
only triplets of MUBs turned up. 

The 12th roots list we arrived at is 

F(0,0) admits 4 candidates for (MUB) 2 . 
F(l/6,0) admits 1 candidate for (MUB) 2 . 
F T (l/6,0) admits 1 candidate for (MUB) 2 . 

This is all. The 4 bases that are MUB with respect to F(0, 0) will be discussed 
in the next section, the basis that is MUB with respect to F T (l/6,0) is an 
enphased version of F(0,0), while the basis that is MUB with respect to 
F(l/6,0) is an enphased version of itself. 
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Figure 2: Here we show all MUB pairs that can be extended to triplets of MUBs, 
and how many triplets there are for each such pair, under the restriction that 
only 24th roots are used. The background for the picture is the set of all known 
equivalence classes of N = 6 Hadamard matrices. 

We redid the entire calculation using 24th roots. This resulted in the 
following additional entries in the list: 

F(l/6, 1/12) admits 4 candidates for (MUB) 2 . 
D(l/8) admits 4 candidates for (MUB) 2 . 

Again this is all; these results are illustrated in fig. 2. In the picture the 
asymmetry between F(xi,X2) and F T (xi,X2) looks quite odd. From section 
3, and from eq. (27) in particular, we know that the unordered MUB pairs 
{l,F(l/6, 1/12)} and {l,F T (l/6, 1/12)} are equivalent, so they should have 
the same number of (MUB) 2 candidates. In section 7 it will be seen that, in 
fact, they do — although in one case none of the candidates is built from 24th 
roots only. 

We redid the calculation using 48th, 72th, and 60th roots; 48 = 2 4 • 3 and 
72 = 2 3 • 3 2 seemed like reasonable things to try, while 60 = 2 2 • 3 • 5 was a 
wild shot. Anyway the list did not increase. It is tempting to think of vectors 
composed of 24th roots as forming a grid in a relevant part of Hilbert space. 
Through a study of the distances between the candidate (MUB) 2 we could 
have refined the grid in interesting places, but we did not pursue this — the 
grid is not very dense. 
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6. MUB triplets including the Fourier matrix 

Let us consider the case when (MUB)! is the Fourier matrix. Here we do not 
have to rely on our own calculations, because complete results are available: 
using the symbolic manipulation program MAGMA, Grassl [14] has com- 
puted all vectors unbiased with respect to both the standard and the Fourier 
basis. The same calculation was in fact done earlier by Bjorck and coworkers 
[15, 16]. There are 48 such vectors altogether, or 24 vectors and their com- 
plex conjugates. Each vector shows up in two different bases, so they form 
16 bases altogether. Using the classification of Hadamard matrices, the 16 
(MUB) 2 turn out to form four groups: 

(i) Two circulant Fourier matrices enphased with 12th roots of unity; "wF' . 

(ii) Two matrices equivalent to F T (l/6,0), and enphased with 12th roots of 
unity; "c;F T ". 

(iii) Six circulant Bjorck matrices C (three with d and three with d), enphased 
with 12th roots of unity; "cjC" . 

(iv) Six Fourier matrices enphased with products of Bjorck's magical number 
d (and d) and 12th roots of unity; u du;F" . 

See eq. (10) for the definition of d. We will refer to the bases in groups (i) 
and (ii) as classical, and to the remaining bases as non-classical, for a reason 
that will transpire. The matrices in groups (i) and (iii) are circulant when 
their columns are multiplied by suitable phase factors. They already include 
all vectors unbiased with the Fourier matrix. Both the ti;F T matrices contain 
three columns each from the two uF matrices. Each matrix duF contains 
one vector from each matrix uC, and vice versa. The four groups of bases 
are invariant under complex conjugation. Two of the bases in uC, namely 
those given by Bjorck's matrix with d and with d, are unchanged by complex 
conjugation. 

The distances between the 16 bases are as follows: The classical bases 
form a perfect square, with edge lengths squared D 2 C = 0.4 (and the enphased 
Fouriers in opposite corners). The distance between any classical and any 
non-classical basis is always given by D\ = 0.92. The distance between any 
member of group (iii) to any member of group (iv) is always D\ = 0.74. 
Finally the distance relations within the two groups of 6 are identical; each 
group contains two equilateral triangles with edge lengths D\ = 0.93. These 
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Figure 3: All Hadamard bases MUB with respect to F (a single point in fig. 2). 
The square consists of (enphased, affine) Fouriers, one of the two David's stars 
of Bjorck matrices enphased with 12th roots, and the other of Fouriers enphased 
with the number d. The numbers refer to the chordal distance squared, D 2 . 

triangles are exchanged if d and d are exchanged. In addition, within each 
group each basis has two bases at D 2 = 0.86 and one at D 2 = 0.88. 

The Bjorck/Grassl results are summarised in fig. 3. There does not exist 
any set of four MUBs including 1 and F, since the distance squared between 
any possible pair of (MUB) 2 never reaches 1. The best we can do with 7 
bases is D 2 > 0.86. This can be compared to random searches; in a sample 
of 20 million bases, chosen according to the measure described in section 4, 
the best results were D 2 > 0.91 for 4 bases, and D 2 > 0.86 for 7 bases. 

It is striking that every vector vector unbiased with respect to 1 and F can 
be collected into some (MUB) 2 . The "doubling" of the non-classical bases 
is striking too. To explain these features we begin by proving that a vector 
is unbiased with respect to both 11 and F if and only if it is a member of a 
circulant Hadamard matrix. This requires us to recount some well known 
facts about the discrete Fourier transform. Given a sequence of complex 
numbers z a , < a < N — 1, we define its Fourier transform to be 
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I N-l 



N b=0 



E <A& ■ (so) 



In condensed notation 



z = Fz z = ¥ ] z . (51) 

If we are given a sequence of complex numbers z a , the column vector whose 
components are z a j \/N is unbiased with respect to the Fourier basis if and 
only if the sequence z a is unimodular, |^ a | 2 = 1, and it is unbiased with 
respect to the standard basis if and only if z a is unimodular. Hence vectors 
that are unbiased with respect to both the standard basis and the Fourier 
basis are in one-to-one correspondence to sequences obeying 

\z a \ 2 = \z a \ 2 = 1 (52) 

for all values of a. Such sequences are called biunimodular. (Standard con- 
ventions for the Fourier transform sit a little uneasily with our conventions 
for MUBs. We decided to live with this, which is why our vectors are made 
from z a rather than z a itself.) 

Biunimodular sequences have an interesting property, that emerges when 
one studies the autocorrelation function 

Y N-l 

7> = T7 2 ~ Zala + h ■ ( 53 ) 

iV a=0 

An easy calculation shows that 

1 N-l 

7. = E W V 6 • (54) 

iV a=0 

Hence, if the sequence is biunimodular it obeys 

lb = 8b,o ■ (55) 

Therefore z a and z a +b, with b fixed and non-zero, are orthogonal vectors. 
A matrix is called circulant if it takes the form 
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z Z N ^i 
Z\ Zq 



Zl 
Z2 



(56) 



ZN-l Zn-2 ■ ■ ■ Zq 



The matrix elements are 



C a b — z a -b ■ (57) 

With this definition a circulant matrix is Hadamard (up to normalisation) 
if and only if the sequence z a is biunimodular. It follows that all vectors 
unbiased with respect to both the standard and the Fourier bases can be 
collected into a set of circulant Hadamard matrices whose columns form 
bases that are MUB with respect to the standard and Fourier bases. 

This explains why all the 48 unbiased vectors are included in some basis 
that can be represented by an Hadamard matrix in circulant form. It remains 
to understand the "doubling" of bases found in Grassl's list. For any circulant 
matrix C, with first column z a , it is easy to check that 

F f C = DF* , D = diag(z , *i, . . . , ^-i) • (58) 

Together with the earlier results, this entitles us to say that a vector is 
unbiased with respect to both the standard and the Fourier bases if and only 
if it is a column of a circulant Hadamard matrix, and every such matrix 
represents a basis that is MUB with respect to both 11 and F. Incidentally 
it also implies that the Fourier matrix diagonalises every circulant matrix. 
Now consider an unordered MUB triplet {1, F, H}. From section 3 we know 
that 

{1, F, H} « {F t , 1, F^H} « {F, 1, F^H} . (59) 

In the last step we multiplied F* from the right with P = F 2 , which is a 
permutation matrix. Now let H = C, a circulant matrix. Using eq. (58) 

{1, F, C} w {F, 1, F f C} w {F, 1, DF 1 "} . (60) 

This is how the second "David's star" of phased Fouriers, in fig. 3, arises. 
The transformation taking one David's star to another also exchanges the 
two enphased Fourier MUBs, while it leaves the pair in group (ii) invariant. 
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For any dimension the conclusions are, first, that every vector unbiased 
with respect to both 11 and F is given by a biunimodular sequence, and 
conversely. Second, such a vector exists if and only if it is a member of a 
circulant basis that forms a MUB triplet with 11 and F. Third, if this circulant 
basis is inequivalent (as an Hadamard matrix) to the Fourier matrix then 
there exists an equivalent MUB triplet whose third member is an enphased 
Fourier. We have no explanation for the two MUB triplets involving members 
of the F T family — from our point of view they arise by accident. 

Counting the number of unbiased vectors, or equivalently biunimodular 
sequences, remains a dimension dependent matter. By means of symbolic 
manipulation programs all biunimodular sequences with N < 8 entries have 
been listed by Bjorck and coworkers [15, 16] (see also Haagerup [9]). For 
N = 6 there are 48 biunimodular sequences (up to multiplication with a 
complex phase), or equivalently 48 vectors that are MUB with respect to 
both the standard and the Fourier bases. This includes 12 classical ones built 
from 12th roots — they were known to Gauss — and 36 non-classical ones built 
from Bjorck's magical number d. 

7. MUB triplets with a block structure 

In section 5 we observed that there were several MUB triplets built from 
24th roots. The ones that cannot be reduced to 12th roots are as follows: 

• Starting from the MUB pair (1, D(l/8)) we have 4 triplets including 
enphased versions of F(l/6, 1/12). They form a semi-regular simplex 
where 4 edges have D 2 = 0.89 and 2 edges have D 2 = 0.95. 

• Starting from the MUB pair (1, F(l/6, 1/12)) we have 4 triplets includ- 
ing enphased versions of D(l/8). They form a regular simplex where 
all the edges have D 2 = 0.93. 

The value D 2 = 0.95 is the largest we have found in our searches for approx- 
imate MUB quartets. 

These MUB triplets have an interesting structure in common. The affine 
family including the Fourier matrix, given in eq. (4), consists of Hadamard 
matrices that are equivalent under permutations to 
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F(x 1 ,x 2 ) « F D 



F 3 
F 3 D 



F 3 

-F 3 D 



D = diag(l,zi,z 2 ) . (61) 



Here F 3 is the 3 by 3 Fourier matrix. This is a kind of twisted version [9, 10] 
of the tensor product matrix F2<S>F 3 .The Di^a family of matrices also has an 
interesting block structure. It was noted by Zauner that they can always be 
written on block circulant form, meaning that they can be brought to a form 
where the matrix consists of four blocks, each of which is a 3 by 3 circulant 
matrix [12]. This turns out to be relevant here. 

Consider the MUB pair (1, F(l/6, 1/12)). Write it as 



(1, F D ) , D = diag(l, 1, cv 6 ) , u = e 2m/M . (62) 

According to our list this pair can be extended to four triplets including an 
Hadamard matrix equivalent to D(l/8). With the members of the pair fixed 
according to the above, the four possible choices for (MUB) 2 are 
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where the Cj are the 3 by 3 circulant matrices 
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In this sense the block structures of the Fourier and Dita families are related. 
We now come to the following little puzzle. We know that 



{1, F(l/6,1/12), D(l/8)}« 
» {[F(l/6, 1/12)]+, 1, [F(l/6, l/12)]tD(l/8)} » (65) 
{1, F T (1/6,1/12), L>P[F(l/6,l/12)]tD(l/8)}, 
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where DP is the product of a diagonal unitary and a permutation matrix. 
Why then does F T (l/6, 1/12) not appear in fig. 2? The answer is that the 
matrix [F(l/6, l/12)] t D(l/8) is not built from 24th roots only. 
What is needed to resolve the puzzle is the number 

h = e 2mc \ cos(27r Cl ) = y|. (66) 

We modified the "24th roots program" to allow also b\ and b\ as entries in 
the vectors. We found the following new MUB triplets: 

• Starting from the MUB pair (1, F T (l/6, 1/12)) we have 4 triplets in- 
cluding enphased versions of F T (ci, 0). 

• Starting from the MUB pair (1, F T (ci,0)) we have 2 triplets including 
enphased versions of F T (1/6, 1/12). 

• Starting from the MUB pair F(ci,0) we have 2 triplets including en- 
phased versions of D(— 1/8). 

The first item in this list contains the "missing" triplets whose existence we 
already knew about — the asymmetries in fig. 2 were due to the 24th roots 
restriction only. 

We have made computer searches to see if we can find other phase factors 
increasing the number of MUB triplets. One such number turned up, viz. 

b 2 = e 2niC2 , tan (2ttc 2 ) = -2 . (67) 

Again we modified the "24th roots program" to allow also b 2 and b 2 as entries 
in the vectors. Then we found the following new MUB triplets: 

• Starting from the MUB pair (1, D(0)) we have 2 triplets including 
enphased versions of F(9/24 + c 2 ,0). The distance squared between 
the latter is D 2 = 0.77. 

• Starting from the MUB pair (1, F(9/24 + C2,0)) we have 120 vectors 
unbiased to both, of which 60 vectors form 10 bases that are enphased 
versions of D(0). They form a polytope where 9 edges end at each 
corner, 6 with D 2 C = 0.93 and 3 with D 2 = 0.78. 
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So we have 6 bases with D 2 C > 0.93. 

An explicit example of a MUB triplet including 62 is 

(l,F D ,D bc ) . (68) 

Here D = diag(u; 9 &2, 1,1), F D was defined in eq. (61), and D bc is a block 
circulant matrix equivalent to D(0), namely 
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(70) 



Thus it appears that many of the MUB triplets in this section can be brought 
to the (twisted product Fourier) — (block circulant Di^a) form. But we did 
not check them all, nor do we have any complete results analogous to those 
available for the (1, F) pair. 



8. What is the most general Hadamard matrix? 

We now take up the story of Hadamard matrices again. In section 2 we 
gave a list of all known Hadamard matrices of order 6, or equivalently of 
all ordered MUB pairs up to equivalences. The question is whether this 
list is complete. To investigate this question one can try a perturbative 
approach to the unitarity equation: multiply the non-trivial matrix elements 
of a dephased Hadamard matrix with arbitrary phase factors, and solve the 
unitarity equation to first order in the phases. The number of free parameters 
that are left when this has been done is an integer known as the defect of 
the Hadamard matrix. It provides an upper bound on the dimensionality of 
any analytic set of Hadamard matrices [11]. 

The defect for the known Hadamard matrices has been computed. For 
the Fourier, Bjorck, and Di^a matrices it equals 4, while it equals for 
Tao's matrix. Closer inspection shows that the defect is likely to be constant 
within each affine family, with possible exceptions at isolated values of the 
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affine parameters [11]. We have also computed the defect for one thousand 
randomly chosen points along the non-afhne family B(#); it was always equal 
to 4. These results show that Tao's matrix is an isolated point in the set of 
all Hadamard matrices. For the remaining matrices the answer is less clear. 
It appears that their defect is always equal to 4, but the defect provides us 
with an upper bound on the dimensionality only. It is not known whether 
the bound is attained, nor indeed whether additional Hadamard matrices, 
not connected to any of the above, exist. 

Let us be fully explicit about the calculation of the defect for Dita's 
matrix D. We focus on this particular Hadamard matrix because its special 
form simplifies the resulting calculations. Starting from the known form of 
the matrix — see eq. (11), with z = 1 — we multiply all non-trivial matrix 
elements with arbitrary phases according to 

D ab -> D ab e ix ° b , 1 < a, b < 5 . (71) 

The unitarity conditions give 15 complex equations, or 30 real equations; 
more than enough to determine the 25 parameters x a b- 

We will try to solve the unitarity equations order by order in the phases. 
If we expand to first order in the phases we obtain the matrix 
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• (72) 



Let us introduce the (unusual) notation x aa , £( a b), £[ a b] f° r the diagonal, 
symmetric off-diagonal, and anti-symmetric parts of x a b. One finds that the 
30 equations split naturally into 

5 equations for x aa 
10 equations for x^) 

15 equations that split into 6 for £[ ab ] and 9 for x^ ab y 
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To first order all the equations are linear. There will be 4 undetermined 
parameters in X[ ab j. One also finds 



x aa = x (ab) = . (73) 

Since the 19 equations for X( ab ) are linear there are no consistency problems. 
It is however worth observing that the general solution of the final 9 equations 
for x {ab) is 

X(ab) = x , (74) 

with x remaining as a free parameter, to be set to zero by the other 10 
equations for X( ab y 

A suitable choice for the free parameters in X[ ab j is £[12], £[13], £[24], £[34]- 
Then one finds 

^[14] = ^[12] + £[34] 
£[15] = £[13] + £[34] 

£[23] = £[13] + £[24] /^x 
£[25] = -£[12] + £[13] 
£[35] = -£[24] + £[34] 
£[45] = —£[12] — £[24] • 

In conclusion, to first order there are 4 free phases, or equivalently the defect 
equals 4. 

We have carried the calculation to third order. However, the Ansatz in 
eq. (71) must then be modified. Set 

71 , _s. D , Ji x ab+yab+Zab + ---) _ 

(76) 

= D ab ^1 - -x 2 ab - x ab y ab + i(x ab + y ab + z ab - -x 3 ab ) + . . )j . 

Here it is understood that y ab and z ab are quadratic and cubic functions of x a b, 
respectively. For obvious reasons we do not write out the resulting matrix 
here. There is an amount of ambiguity in the definition of the higher order 
terms. This we resolve by imposing 

V[12] = V[13] = V[24] = V[34] = , (77) 



27 



and analogously at all higher orders. This can always be achieved through a 
redefinition of the free parameters in x a b- Once this is done the expansion in 
the exponent of eq. (76) is unambiguous. 

At each order higher than the first, there will be 30 linear inhomogeneous 
equations for the 25 — 4 new parameters that we introduce, with the inhomo- 
geneous terms made up from products of known lower order contributions. 
Given eq. (77), the anti-symmetric and diagonal parts will be completely 
determined. There remain 10 + 9 equations for the symmetric off-diagonal 
part, and it is not a priori clear that a solution exists. 

At second order one finds 

V[ab] = (78) 

2|/ll = ~ x \\ X 31 X W ~ X 51 
2?/22 = ~ x 12 ~ X 32 X 42 X 52 

2^/33 = xj 3 - x 2 23 - x 2 43 + x 2 3 (79) 

2|/44 = 2^14 4" 3^24 "^34 X 54 
2?/55 — ~ x 15 + x 25 + ^35 — X 45 • 

The last 9 equations imply 

V(ab) = V , (80) 

for some y. Detailed examination shows that the remaining 10 equations are 
consistent with this and with each other. The solution can be written as 



Ay = 2y n + 2y 22 - (x 13 - x 23 ) 2 + (xu - x 24 ) 2 - (x 15 - x 25 ) 2 ■ (81) 

The expansion to second order is therefore consistent. 
At third order one finds 

Zaa = , (82) 

together with rather complicated expressions for the part of left unde- 
termined by the condition analogous to (77). We do not give them here. The 
first 10 equations for the symmetric off-diagonal part is solved by 

Z { ab) = . (83) 
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Thus z a b is, like x a b, anti-symmetric in its indices. We are then left with 9 
consistency conditions on this solution. Again detailed examination shows 
that they are obeyed, hence the expansion is consistent to third order. 

In conclusion we have proved that Diva's matrix admits an expansion 
to third order in four arbitrary phases. Unfortunately, to second and third 
order this result depends on cancellations that we do not understand — we 
have simply checked that they do occur. Nevertheless it seems to us likely 
that the expansion exists to all orders. This suggests that the parameter 
space of allowed Hadamard matrices has a four dimensional component. 

Some solutions are known to all orders. First of all, if the expression 
D a b D ao e tXab is exact, that is if all higher order terms in the exponent 
vanish, we have an affine family. It is known that Dita's matrix admits 5 
permutation equivalent affine families [11]. They result from the five choices 



i: 
ii: 
iii: 
iv: 

v: 



(Z[12],Z[13],£[24],Z[34]) = (X, 0,0,0) 

(x [12 },x [13] ,x m ,x [34] ) = (0,2;, 0,0) 

(X[12] , X[13] , X [2 4] , X[ U ] ) = (0, 0, X, 0) (84) 

{X[l2],X[ l3 ],X[24],X[ 3A] ) = (0,0,0,2;) 

(X[12] , X[i 3 ] , X[24] , X[ 34 ] ) = (X, X, -X, ~x) . 

In fact our choice of independent variables was made with this result in view. 
In these cases we have y a b = z a b = 0, and similarly for all higher order terms, 
so these solutions are somewhat trivial. Fortunately at least one non-affine 
family, with non-vanishing higher order terms, is also known to all orders. 
This is the one-parameter family B(0), which passes through representatives 
of the Dfja equivalence class no less than three times. A solution to our 
equations, equivalent to B(0) to third order, is 

(Z[12],Z[13],Z[24],Z[34]) = (-X,-X,2x,x) . (85) 

The higher order terms in the exponent of eq. (76) are now non-vanishing. At 
the same time we observe that, even if it does exist to all orders, this solution 
must encounter convergence problems in order to reproduce the limits on the 
phase 9 in the Beauchamp-Nicoara family. Clearly such a problem can occur 
in the exponent of eq. (76). 

Our calculation, together with several different one-parameter solutions 
known to all orders (one of which is a non-affine solution), does support the 
conjecture that the set of Hadamard matrices is four dimensional in a neigh- 
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bourhood of the Dfya matrix. We understand that Beauchamp and Nicoara 
have numerical evidence for the same state of affairs in a neighbourhood of 
the Bjorck matrix [27]. At the same time the existence of Tao's matrix shows 
that the set of Hadamard matrices must be disconnected. The Bjorck and 
Di^a matrices are connected to each other, but we have no evidence to show 
how the Fourier matrix fits into the picture. 

9. Summary 

Mutually unbiased bases live in the landscape of complex Hadamard matri- 
ces. A complete classification of the equivalence classes of N = 6 complex 
Hadamard matrices, or equivalently of the set of all MUB pairs, is however 
not available. If we insert arbitrary phases in the known examples, we find 
that four parameters are left undetermined to first order. The one exception 
is Tao's matrix. For Diva's matrix we solved the unitarity equations to third 
order in four arbitrary phases. This does not constitute a proof that the pa- 
rameter space has a four dimensional component, but we may perhaps quote 
Lewis Carroll at this point: "I have told you thrice. What I tell you three 
times is true." 

A set of MUBs can always be described as the standard basis together 
with a set of Hadamard matrices. We searched for such sets, mostly under 
the restriction that all vectors be composed from rational roots of unity, 
but we also gave three examples of non-rational and "MUB-friendly" phase 
factors. The restriction to rational roots is severe, but the advantage is that 
systematic searches can be made. 

We found only triplets of MUBs. A rich source of triplets is obtained 
by choosing the Fourier basis as one of the members; the resulting structure 
can be largely understood through the discrete Fourier transform. Other 
triplets, with a special block structure, were also found. To add structure 
to the problem we introduced a natural distance between bases in Hilbert 
space. It attains its maximum value when the bases are MUB. Whenever a 
pair of MUBs (including the standard basis) could be extended to a triplet 
in more than one way, we computed the distance between the possible third 
members. In this way we found sets of 4 bases where all distances but 
one are maximal. The highest figure we obtained for the last distance was 
D\ = 0.95, significantly closer to the maximum value 1 than is the average 
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distance squared between bases, {D%) = 6/7. 

No-go statements for complete sets of 7 MUBs in N = 6 have been made. 
They tend to exclude sets with a high degree of symmetry [13] or with oth- 
erwise special group theoretical properties [28], as well as sets constructed 
using generalisations of methods that work in prime power dimensions [29]. 
We are unsurprised by our failure to find any complete set — but if our conjec- 
ture concerning the set of all Hadamard matrices is true, the no-go theorems 
are unlikely to apply in general. 

So what is the best approximation to a complete set of MUBs that one 
can find in 6 dimensions? Indeed, what is the maximal number of MUBs one 
can have in 6 dimensions? The jury is still out. 

Note added in proof: After submission of this paper Matolcsi and Szollosi 
[30] found a new family of 6 6 Hadamard matrices. And in Ref. [31], Tadej, 
Zyczkowski, and Slomczynski study the defect of unitaries, not necessarily 
Hadamard. 
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